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I. INTRODUCTION 


The existence of neutral non-uniform black string solu¬ 
tions was originally postulated by Gregory and Laflamme 
when they discovered the perturbative classical instabil¬ 
ity of the uniform strings ||, 0. By compactifying 

the axis of rotational symmetry, the solutions can be sta¬ 
bilised above a critical mass. At the critical mass, the 
marginal mode is static indicating a new branch of solu¬ 
tions with spontaneously broken translational invariance 
along the compact direction. The end state of the decay 
for the unstable compactified uniform strings was sup¬ 
posed to be a compact black hole solution, with spherical 
spatial horizon topology. However, Horowitz and Maeda 
Ig, ^ argued that the horizon can not ‘pinch off’ in a fi¬ 
nite affine time, and concluded that the end state should 
instead be a non-uniform string. The non-uniform so¬ 
lutions connected to the Gregory-Laflamme (G-L) criti¬ 
cal point were then constructed in third order perturba¬ 
tion theory by Gubser 0 , who found that their mass in¬ 
creased with non-uniformity, whilst the asymptotic com- 
pactification radius is held constant. In a recent work 
1^, ^ we have shown that the mass of these non-uniform 
strings is in fact always greater than that of the criti¬ 
cal string. This implies they can not be the product of 
decay of an unstable uniform string, and the mysterious 
dynamics of the decay remains an open question. Nu¬ 
merical work in progress may shed light on the product, 
although the evolution time is not long enough at present 
1^, 1^ . Other interesting issues are links between classi¬ 
cal and thermod yna mic stability for the uniform strings 
0 0 ^ 0 , the classical instability in related 

spacetimes |i6| , 0 0 0, and the phenomenological 
implications these strings may have in theories with com¬ 
pact extra dimensions ©Ell- 

Some analytic results are known for non-uniform 
strings |^, 0. However, the Weyl ansatz which solves 
4-dimensional static axisymmetric gravity does not gen¬ 
eralise to higher dimensions to give spherical symmetry 
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. Interesting work by Har- 
has conjectured an ansatz for the 
general string/black-hole solution that reduces the num¬ 
ber of unknown metric functions from three to one, al¬ 
though they could not show its consistency. Using this 
ansatz they conjectured a new branch of non-uniform 
string solutions, unrelated to those connected to the G-L 
critical point. We term these ‘conventional’ non-uniform 
solutions G-L strings, and it is these that mostly concern 
us here. We term the solutions conjectured to exist by 
Harmark and Obers the H-0 strings. 

Using Morse theory, Kol [© has conjectured that 
the G-L string branch of solutions should connect to 
the black hole branch (assuming it actually exists - see 


31 ) via a spatial horizon topology changing solu¬ 


tion. In this Letter, we use our numerical solutions for 
highly non-uniform G-L strings to study the geometry 
of the horizon, as a function of the minimal sphere ra¬ 
dius, whilst keeping the asymptotic radius of compact- 
ification fixed. Following Gubser we will use the quan¬ 
tity, A = i where TZmax is the maximum 


radius of a 3-sphere on the horizon and TZmin is the min¬ 
imum, being the radius at the ‘waist’. Thus A is zero for 
the uniform black string. We briefly review the numer¬ 
ical method, and then provide evidence that as A —> oo 
these solutions are consistent with Kol’s conjecture, that 
A = oo is the horizon topology changing solution, which 
connects with the black hole branch. Finally we demon¬ 
strate that our solutions can always be written in the 
form of the Harmark and Obers ansatz, thereby showing 
consistency of this ansatz, at least for G-L strings. We 
comment on the implication of this for the existance of 
H-0 strings. 


II. CONSTRUCTING SOLUTIONS 


The problem of constructing compactified non-uniform 
string solutions is an elliptic one. We wish to have an 
asymptotic product geometry of flat space and an S^. 
Near the axis of rotational symmetry we wish to im¬ 
pose horizon boundary conditions, with a fixed amount 
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of ‘wiggliness’. In addition we wish the solutions to be 
periodic in the axis direction. The scale invariance of the 
vacuum Einstein equations means that finding the solu¬ 
tions for a fixed size allows all other solutions to be 
generated simply by scaling. 

We first developed numerical methods in to fully 
solve the geometry of a star near its upper mass limit on 
a Randall-Sundrum brane [^, Later in we ex¬ 
tended these methods to the case of the black string. For 
technical reasons, we consider the 6 dimensional black 
string solution, rather than the 5 dimensional one exam¬ 
ined by Gubser. As shown in |^, these have the same 
qualitative thermodynamic behaviour as 5 dimensional 
strings, to leading order in A. We take the metric, 

=- - — + dz'^) + {m + r'^)dnl 

m + r‘^ 

( 1 ) 

where A, B, C are functions of r, z and 2 : = [0, L] is an 
interval as the string is wrapping an S^. This form can 
always be locally chosen. When A,B,C = Q this reduces 
to a uniform black string. As found in [H, the diag¬ 
onal form, and residual invariance under r, z conformal 
transformations results in several important simplifica¬ 
tions. Firstly, the ‘interior’ Einstein equations, 
and (Cj. -I- take the form, 

= srcx, (2) 

for Xi = {A, i?,G}, where = 9^ -I- 9^, and the 
sources srcxj depend non-linearly on all the Xj, drXj 
and dzXj. The form of these equations, or equivalently 
the Lagrangian from which they are derived, allows stan¬ 
dard elliptic numerical methods to be employed to solve 
for A, B, G given elliptic boundary data. Note that the 
superspace metric has Lorentzian signature, so this is 
an extremisation problem, rather than a minimisation 
one. The second crucial observation is that the remain¬ 
ing Einstein equations, G''^ and (G’'^ — G^^,), weighted by 
the measure g = ^ — det obey an elegant Cauchy- 
Riemann relation, 

d,{gG\) + dz[\{G\-G\)) =0 

dz{gG\)-dr[^{G\-G\)) =0 (3) 

provided that the interior equations are satisfied. If, for 
example, we satisfy gG"^^ on all the boundaries, we need 
only satisfy g{G''^ — G^^) at one point and both con¬ 
straints must then be zero everywhere in the interior. 
Asymptotically we choose A,B,C 0, and thus L is 
the radius of the compactification. We also enforce the 
{G'^j. — G^z:) constraint there, and the G'^^: constraint is au¬ 
tomatically true as 2 : dependence decays exponentially at 
large r with these boundary conditions. The third impor¬ 
tant feature resulting from the r, 2 : conformal invariance 
is that we are free to choose the boundaries, where we fix 
data, to be at any coordinate location. Thus we need not 
include additional information to encode boundary posi¬ 
tions. We take periodic boundary conditions at z = 0, L, 


and choosing the horizon to be at r = 0, compatible with 
our choice of background string solution, implies that 
A,B,C will be finite at the horizon. Horizon boundary 
conditions are imposed by regularity We find [|| that 
we must impose A^r = 0, G ^ = 0 p6| and the constant 
temperature condition, 

(dzA - dzB) \r^o= 0 (4) 

which ensures G’’^, is satisfied, and must be integrated 
along the horizon, giving rise to one parameter, B^ax at 
z = L, as the integration constant, that determines the 
deformation from the uniform string. From our numeri¬ 
cal solutions, we find a one-to-one relation between A and 
ddmax ■ 

Relaxation of the interior equations, with the bound¬ 
ary conditions described above, is continued until a con¬ 
vergent solution is found. Using elementary numerical 
methods and moderate resolution and computer time, so¬ 
lutions were found for A < 4. In figure ^ we show the spa¬ 
tial horizon geometry of the most non-uniform solution, 
embedded into Consistency checks were performed in 
[^; direct evaluation of the constraints, comparison with 
third order perturbation theory at small A, varying L 
(which gives scale equivalent solutions). The asymptotic 
ADM mass of the solutions is plotted in figure || against 
A. This can be computed from the asymptotic behaviour 
of the metric, or from integration of the First Law, and 
the two give good agreement. For A ~ 4, the difference 
in the values is approximately 5%. For a detailed discus¬ 
sion of errors, see |^. It is difficult to assess systematic 
errors in the method, but it is likely that other data plot¬ 
ted here have similar magnitude errors. From this figure 
H, we concluded that the mass of these non-uniform G-L 
strings is always larger than the critical string, and for 
A —> 00 it tends to approximately twice the value of the 
critical string. 

III. KOL’S CONJECTURE 

We now consider how A becomes large. For our so¬ 
lutions, the maximum value for the metric function G 
occurs at z = 0, and the minimum is at z = L, where 
B = Bjnax- This gives rise to the horizon 3-sphere radii 
Tlmax and TZmin, plotted in figure ^ as a function of 
1/(1 -I- A), for fixed asymptotic radius, and normalised 
by the critical uniform string horizon radius. Clearly for 
A to become large, either TZmin 0, or TZmax 00 , 
or both. It is immediately clear from the figure that 
A becomes large due to TZmin shrinking, apparently as 
1/A for large A. The earlier horizon embedding, figure 
□, graphically illustrates this, the critical uniform string 
with the same asymptotic size having unit horizon 
radius. Thus we find good evidence that TZmin ^ 0 as 
A ^ 00 . 

Taking the black hole branch of solutions, for the same 
asymptotic compactification radius, and increasing their 
mass from zero, the black holes will appear less and 
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FIG. 1: Horizon geometry of non-uniform string found with 
A ~ 4. For the same asymptotic radius, the critical uni¬ 
form string has unit horizon radius. This indicates that the 
large value of A is mainly due to a cycle shrinking, suggesting 
the possibility of a topology change for A —> oo. 





FIG. 2: Plot of mass, normalised by the critical uniform 
string mass, against A for fixed asymptotic compactification 
radius. The data points indicate actual solutions. Three reso¬ 
lutions are shown superposed, the highest resolution allowing 
solutions to be found with A ~ 4. 


less like 6 dimensional Schwarzschild solutions. At some 
point, the branch of solutions might cease, there could 
be a point where the horizon topology changes H), 
the solutions might decompactify along the axis, as oc¬ 
curs for the extremal charged black hole solutions, or the 
black hole could become highly extended only in the ra¬ 
dial direction j2^, If, as Kol suggests, there is a point 
of topology change and the solutions connect to the G-L 
non-uniform branch at A = oo, so TZmin = 0, we would 
expect a finite TZmax- On the other hand, if this is not 
the case, the geometry of the G-L strings may do some¬ 
thing quite different as A ^ oo. An example would be 


if Harmark and Obers are correct, and it is a different 
branch of non-uniform solutions - the H-0 strings - that 
connect to the black hole solutions, rather than the G-L 
solutions. 


Normalised R . and R 
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FIG. 3: Plot of maximum and minimum horizon radius, nor¬ 
malised by the critical uniform string radius, against 1/(1-|-A), 
for fixed asymptotic size. We see Rmin oc 1/A for A —> oo 
and Rmax remains finite. 

We see from figure || that TZmax does appear to tend 
to a constant value as A —> oo. Thus by considering the 
extrema of the horizon, we find consistency with the idea 
that the black hole solutions could connect to the non- 
uniform branch. Let us now consider other geometric 
quantities, to see if pathologies develop that could spoil 
this picture. 

Firstly, whilst we fix the asymptotic length of the S^, 
we do not a priori know the proper length along the hori¬ 
zon, Lhoriz, found by integrating at r = 0. We plot 
this in figure^, and see that it increases, but appears con¬ 
sistent with tending to a constant as TZmin 0, again 
being compatible with the solution tending to the lim¬ 
iting black hole solution which can only just ‘fit’ into 
the circle. Note that another possible option would have 
been finding this horizon length became infinite in the 
limit, which obviously would not have been compatible 
with Kol’s conjecture. 

Thus we have considered coordinate invariant quan¬ 
tities associated with the B, C metric functions. The 
remaining metric function A forms the coordinate invari¬ 
ant horizon temperature, as . As shown in |^, and 
in figure ||, the temperature again appears to tend to a 
constant. This is crucial as we expect this temperature 
to remain finite if there is a continuous interpolation be¬ 
tween black hole and string solutions. 

Note that finding finite values for both the above quan¬ 
tities is encouraging from a numerical point of view, as 
at z = L, where the horizon radius is a minimum, we find 
both e^ and e^ appear to diverge as A —> oo. Obviously 
we expect the coordinate system to become singular at 
the topology changing point, and it is good news that 
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such coordinate invariant quantities remain well behaved 
there. 


Normalised length of horizon 



FIG. 4: Plot of proper length of horizon against Rmin, nor¬ 
malised by the critical uniform string value, for fixed asymp¬ 
totic radius. 


Normalised horizon temperature 
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FIG. 5: Plot of horizon temperature against Rmin, nor¬ 
malised by same the value for the critical uniform string, for 
fixed asymptotic size. 

To summarise, several geometric quantities appear to 
have limits as TZmin 0 compatible with interpolating 
to black hole solutions, as conjectured by Kol. We stress 
that this is by no means a proof. However, the conjec¬ 
ture could easily have been falsified if any of the tested 
quantities had not behaved in a finite way as A —> oo, 
so we do regard it as supporting evidence. Assuming 
Kol’s picure holds, the values of Rmax, Lhorizon and the 
horizon temperature in the above figures, extrapolated to 
A —> oo, can then be taken as testable predictions for the 
same quantities measured in the largest mass black hole 
solution. If, for example, our numerical methods can be 
modified to explore the black hole branch, we may then 


quantitatively test whether the branches could join, even 
if the topology changing solution can not be numerically 
constructed explicitly. 

IV. HORIZON GEOMETRY AS A ^ oo 

We now consider the geometry at the minimum radius 
3-sphere of the horizon. A simple scalar measure of cur¬ 
vature is the Kretschmann invariant, K = RaPfiu- 

Since we are interested in its behaviour at this extremal 
point of the horizon, r = 0, z = L, and dr^zXi = 0 there, 
we may expand to quadratic order in r, z, 

Xi ~ a;(o)i + -|- X(2)i'r^ (5) 

where = {a^), cq)} for Xi = {A,B,C}, and 
z = z + L. Using the interior equations we may solve 
X( 2 )i in terms of a;(o.i)i) and the constant temperature 
condition, resulting from the constraints, yields a(i) = 
6 ( 1 ). We find, 

AT = 48 -|- 

( 6 ) 

at the waist r = 0 , z = 0 , choosing units where m = 1 . 
We now know that for this minimal cycle, C(o) ^ —oo as 
A —> oo (as = Rmin —> 0). From the form of the 
equation above, it is easy to see that K must diverge in 
this limit, independent of the behaviour of 6 (q) and C(i). 
Thus we have shown the curvature at the waist must 
become singular, provided the solutions exist for large 
A with the boundary conditions we impose, so that the 
near waist behaviour (i remains true. Assuming that 
A remains finite away from the horizon then implies the 
A = oo solution has a naked singularity at the waist. This 
again agrees with Kol’s picture, who has also suggested 
that this nakedly singular geometry near the waist may 
be a cone metric, which we aim to test in future work. 

V. THE HARMARK-OBERS ANSATZ 

Finally, we discuss the metric ansatz proposed by Har- 
mark and Obers, applied to 6 -dimensional neutral black 
string/hole geometries, namely, 

= — 2^-2 df+e‘^^{df^+e~^^dz'^)+e^^{rl+f^)dnl 
Tq +r 

( 7 ) 

where M, K are functions of r, z, and tq is a constant. In 
Ip^ , Earmark and Obers constructed this ansatz based 
on a coordinate system that ‘interpolates’ between the 
black hole and uniform string solutions. They also find 
that K can be algebraically eliminated in terms of M. 
This therefore gives an ansatz containing only one free 
function, M, rather than the usual 3 generally required 
by static axisymmetry, as in our metric ( 0 )- 
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Of the five independent Einstein equations discussed 
above, one linear combination vanishes, one is used to 
determine K, leaving 3 remaining equations for M. Har- 
mark and Obers were unable to show that these could 
be solved consistently, but gave suggestive evidence by 
considering the consistency of the ansatz at the horizon, 
and to second order in an asymptotic expansion in f. 

In fact it is easy to show the 3 Einstein equations may 
be locally consistent, due to the 2 independent compo¬ 
nents of the contracted Bianchi identities. However, the 
question of global consistency appears to depend on the 
explicit existence of solutions. In this section we exhibit 
the coordinate transformation that takes our metric co¬ 
ordinate system (|^ into the Harmark-Obers coordinate 
system. Our metric contains sufficiently many metric 
functions to parameterise the general solution locally, 
and thus this coordinate transformation ‘locally’ demon¬ 
strates the consistency of the Harmark-Obers ansatz. 
However, our numerical solutions indicate that our coor¬ 
dinates are good globally for the G-L strings (with finite 
A), and therefore shows the ‘global’ consistency of their 
ansatz for these solutions. 

Taking the coordinate transformation, 

f = f{r,z) (8) 

z = h{r, z) 

then gives the conditions, 

drf =+e-^^d,h (9) 

dj =-e-^^drh 

for the transformation of (j^) to Harmark-Obers form. 
In addition, we require that after the transformation the 
lapse is the specific function in (|^). The task is then to 
understand what /, h solve these equations, and consis¬ 
tently gives only this f dependence in the lapse. For this 
lapse to take the desired form, we would require, 

f2 =_ ^0^ ^ _ (10) 

m + r^(l — 

To check whether this can be a solution to the coordinate 
transformation conditions (^, we may eliminate h from 
these to give, 

V^f + 3idrfdrK + dJd,K)^0 ( 11 ) 

Substituting our required / above into this condition, 
gives a second order equation for A which is exactly the 
interior equation, = src^ of (^. Thus the / re¬ 

quired by the Harmark-Obers ansatz is indeed a consis¬ 
tent solution of this coordinate transform, and thus we 
may transform our metric into the Harmark-Obers form. 

We now consider the positions of the horizon, and pe¬ 
riodic boundaries. In our solutions, we have chosen them 
to be at r = 0 and z = 0,L, respectively. We require 
dzXi = 0 at the periodic boundaries. A^B,C are finite 
on the horizon, and go to zero at large r. Together, these 


facts imply that in the f, z coordinates, the horizon is 
at f = 0, and the periodic boundaries are at constant z. 
With a suitable choice of data for h in the transformation 
above we may pick ^ = 0 to be one periodic boundary. To 
find the second periodic boundary location, say z = L, 
one must then use K and (^ to give, 

dzh =[m + r{m + r'^)drA) (12) 

At large r, A, C only have r power law dependence, the 
z dependence having decayed away exponentially. One 
integrates to find h{r, z) = z at large r, which gives L at 
z = L. As A, H, C ^ 0 as r ^ oo so M, K also vanish 
there and thus L is the proper asymptotic circle length. 
It is now important that we may choose ro so that L = L 
ie. choosing tq correctly ‘undoes’ the global scaling that 
will generically occur in this coordinate transformation. 
What is crucial is that for fixed m,L = L, this means 
that ro is a function of the specific solution A, B, C, and 
therefore of our parameter A. 

Plotting the function / shows its r derivative is every¬ 
where positive, within these boundaries, for the solutions 
available with A < 4. Then (^ similarly shows that the z 
derivative of h will be positive everywhere. This implies 
that for the G-L strings tested, no pathologies of the f, z 
H-0 coordinate system arise in the interior. 

Thus the Harmark-Obers ansatz appears to be entirely 
consistent for the G-L branch of solutions. Indeed taking 
ro as the required function of A, the periodic boundaries 
are then at constant z = 0, L, and the horizon at f = 0, 
which is the same as considered in [^. From inspection 
of the Einstein equations resulting from the interpolating 
H-0 ansatz, Harmark and Obers proposed that the black 
hole solutions went through a topology changing solution 
into the H-0 strings, which they thought to be distinct 
from the G-L strings. They parameterised their solutions 
exactly using rp. Inverting the arguments above, we now 
could think of our G-L strings as being parameterised by 
ro rather than A. 

To summarise; we have shown that A —> cx) appears to 
be a topology changing limit for the G-L strings, have 
shown they can be consistently expressed in the H-0 
ansatz, with the same boundary conditions and bound¬ 
ary locations as discussed by Harmark and Obers, and 
furthermore have shown how their parameterisation of 
solutions in terms of rg fits with the one in terms of A. 
In addition, Harmark and Obers also predict a nakedly 
singular topology changing point, as Kol does, and as 
we give numerical evidence for here. It is then intrigu¬ 
ing to consider whether, rather than being a new distinct 
branch of solutions, these H-0 non-uniform solutions are 
simply the G-L strings. This would certainly appear to 
be the simpler outcome. Then Harmark and Obers ar¬ 
gument for the continuation of the black hole solutions 
into some non-uniform string solution appears to be con¬ 
sistent with Kol’s expectations. 

If this is the case, the Harmark-Obers ansatz might 
provide important clues in how to construct the G-L 
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string solutions through the topology changing point, 
where an interpolating coordinate system must be de¬ 
vised which still allows elliptic numerical methods to be 
used. 


VI. CONCLUSION 

Using our recently developed elliptic numerical meth¬ 
ods, we present evidence supporting Kol’s conjecture, 
namely that the compactified G-L nonuniform string so¬ 
lutions connect to the black hole branch of solutions as 
A —> oo. The evidence is based on examining several geo¬ 
metric quantities which must remain finite as the minimal 
horizon sphere shrinks to zero radius, if the conjecture is 
to be true. If true, we may now make quantitative pre¬ 
dictions about the mass and geometry of the maximal 
black hole solution. We find the curvature at the waist 
appears to diverge in the A —> oo limit, implying a naked 


singularity for the topology changing solution. We have 
also shown consistency of the Harmark-Obers ansatz for 
these G-L string solutions, which we see are contained in 
this ansatz. Following from this, we suggest that there 
are no new distinct nonuniform H-0 strings, as conjec¬ 
tured by Harmark and Obers, and the nonuniform strings 
they consider are simply the G-L strings, tying in their 
work with Kol’s picture. 
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